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Abstract 

Refraction of high frequency noise by mean flow gra- 
dients in a jet is studied using the ray-tracing methods 
of geometrical acoustics. Both the two-dimensional (2D) 
and three-dimensional (3D) formulations are considered. 
In the former case, the mean flow is assumed parallel and 
the governing propagation equations are described by a 
system of four first order ordinary differential equations. 
The 3D formulation, on the other hand, accounts for the 
jet spreading as well as the axial flow development. In this 
case a system of six first order differential equations are 
solved to trace a ray from its source location to an observer 
in the far field. For subsonic jets with a small spreading 
angle both methods lead to similar results outside the 
zone of silence. However, with increasing jet speed the 
two prediction models diverge to the point where the par- 
allel flow assumption is no longer justified. The Doppler 
factor of supersonic jets as influenced by the refraction 
effects is discussed and compared with the conventional 
modified Doppler factor. 

1. Introduction 

Propagation and refraction of sound emitted from 
convecting multipole sources embedded in a shear flow has 
its application in jet noise studies. The shrouding effect 
of the mean flow and its contribution to the directivity of 
jet noise has been studied by numerous researchers in the 
past. It is well known that the operator part of Lilley’s 1 
equation contains the effect of refraction as well as the 
source convection, with the right hand side of the equa- 
tion representing, at least to a large degree, the source 
term only. For axisymmetric jets, this equation can con- 
veniently be written in cylindrical coordinates and with 
appropriate transformations a turning point problem is 
formulated. 
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Thus, Mani 2 ’ 3 studied the mean flow interaction of round 
jets at subsonic eddy convection velocities for plug-flow 
jets. The source term in his model is a quadrupole con- 
vecting along the center-axis of the jet. He further ar- 
gues that the field of sources located symmetrically about 
the jet axis is well represented by that of one source on 
the axis. With these assumptions Mani avoided consid- 
erable mathematical difficulty and obtained good agree- 
ment with data for both cold and hot jets except at high 
frequencies. 

The high frequency solution is believed to provide a 
reasonably good approximation to the refraction problem 
even at Helmholtz numbers as small as one 4,5,6 , and is of 
particular interest in the prediction of the noise directiv- 
ity of high speed jets. A high frequency asymptotic so- 
lution for arbitrarily located convected singularities in an 
axisymmetric parallel flow was derived by Balsa 4 . In his 
model, the turbulent properties of the jet and its acoustic 
radiation are coupled through a directivity factor for each 
of the quadrupoles contained within the turbulent eddy. 

To study the radiation pattern of non- axisymmetric 
jets, Goldstein 7 derives a circumferential directivity factor 
in the high frequency limit. For aninviscid parallel flow, 
he obtains a Green’s function for a convecting monopole 
source when the Mach number and temperature are ar- 
bitrary functions of the cross-flow coordinates. The rays 
of geometrical acoustics (GA) are traced by following the 
projection of each ray on a cross-stream plane. Following 
Goldstein’s work, a high frequency Green’s function for 
convecting multipole sources in a spreading jet was devel- 
oped by Durbin 8 and used to predict the noise directivity 
of quadrupole sources convecting along the centerline of 
an axisymmetric jet. It is generally accepted that the high 
frequency acoustic sources he in the neighbourhood of the 
shear layer development, which for an axisymmetric jet is 
simply a ring about the jet axis. However, by studying 
an on-axis source point, Durbin 9 was able to simplify the 
directivity factor significantly and show some important 
features of the refraction problem. 

Implicit in the use of the high frequency approxima- 
tion is the assumption that the acoustic wavelength is 
shorter than the characteristic length of the mean flow. 
In addition, if turbulence length scale is small relative to 
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the wavelength at which noise is produced, Lilley’s equa- 
tion suggests that the sound field is the result of radia- 
tion by independent correlation volume elements within 
the jet, each moving with the local convection velocity. 
A convecting ring source problem was done in 10 by first 
considering an off-axis point source, leading to a 3D ray- 
acoustic problem, and then carrying out a circumferen- 
tial integration for an axisymmetric jet. In addition, the 
mean flow velocity profile was modeled to account for the 
streamwise development of the jet from a top-hat, near 
the exit, to a Gaussian, representative of the downstream 
fully developed flow. The size of zone of silence was found 
to be sensitive to jet spreading angle as well as the mean 
flow gradients surrounding any source point. Further, it 
was shown that the Doppler factor is influenced by the 
refraction effects and varies from point to point on a ring 
source. 

The purpose of the present paper is to further our un- 
derstanding of the refraction of high frequency jet noise 
and its application to jets numerically computed from 
CFD codes. To this end, both 2D and 3D ray-acoustic for- 
mulations are considered and compared. The governing 
equations are briefly discussed in §2 with details relegated 
to the references. Section 3 describes two numerical ex- 
amples. In the first example, the jet is modeled in closed 
form and the directivity angles and the sound pressure 
level are worked out using both 2D and 3D techniques. In 
order to extend the problem to supersonic jets, a second 
example is proposed in which the mean flow is calculated 
numerically. The jet under consideration is issued from 
a convergent-divergent nozzle at its design condition with 
a jet exit Mach number of Uj/a^ = 2.16. We further 
require a shock-free flow at this point, as the propagation 
of rays across shocks discontinuities will require special 
treatment that is beyond the scope of the present work. 
The difficulties arising from the computation of the direc- 
tivity factor of supersonic jets near the boundary of zone 
of silence are discussed. The Doppler factor in the di- 
rection of eddy Mach wave emission is compared with the 
conventional modified Doppler factor based on the line-of- 
sight method 11,12 . Finally the spherical directivity of an 
off-axis isotropic quadrupole convecting at various loca- 
tions along the jet is calculated and compared with data. 
This is followed by summary and conclusions in the final 
section. 


2. Governing Equations 


radiated by a typical correlation volume element. This 
is essentially done by Lighthill’s equation. However, the 
interaction of the mean flow and its effect on the direc- 
tivity pattern of the jet noise can be accounted for in the 
context of Lilley’s equation; where the mean flow gradi- 
ents as well as the source convection are moved to the 
operator part of the equation. In the high frequency ap- 
proximation, as in geometrical optics, sound propagates 
along rays. Hence, the radiation pattern for each convect- 
ing multipole source is the result of the variation of the 
ray-tube area surrounding that ray. Therefore, rays emit- 
ted from each sound producing region of the flow need to 
be traced to an arbitrary observer point. 

2*1 Two-Dimensional Geometrical Acoustics 


For an inviscid parallel flow whose velocity U and 
sound speed a are arbitrary functions of the cross- 
stream coordinates X* , a pair of second order differential 
equations 7 are solved to trace the rays of GA 


sO = v. ? (x,, 

q 2 = $ 2 (X t ) - cos 2 Oca 


*(X t ) = 


1 


C(Xt) 


-M(X,)cose„ 


M(X t ) = 


U(X t ) 


a(X t ) ’ 

subject to the initial conditions 


C(X t ) = 


«(X«) 


(1) 


(X«) = (X t )„, = (cos 6, sin 6) at s = 0. 


In equations (1), s is the arc length along the projection 
of the ray on a cross-stream plane (figure 1), 6 is the radi- 
ation angle at the source, (9 , <p) are polar and azimuthal 
angles in spherical coordinates and subscripts o and oo re- 
fer to source and observer locations respectively. It is pos- 
sible to solve (1) numerically as a system of four first order 
differential equations and continue along each ray until it 
becomes a straight line where <poo is determined as a func- 
tion of emission angle <5. Thus for any pre-selected value 
of polar observation angle 0«> , the spreading or focusing 
of acoustic rays is expressed as the azimuthal directivity 
factor d6/d<j>oo' The far field mean square pressure pro- 
duced by point multipole sources convecting downstream 
at convection Mach Me = Ue/^oo is given as 7 


It is well known that in a fine-grained turbulence, the 
retarded time across a turbulent eddy volume element is 
negligible. At subsonic convection Mach numbers, the ra- 
diation field arriving from each source volume element is 
independent of that due to any other region of the flow 
and the noise field is predicted by estimating the sound 
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Vi = q 0 cos 6 , v 2 = 6, ^3 = cos 6^ 


where is the radiation field due to a monopole 
source of unit strength, P? >1 and Pq.. denote the fields due 
to dipole and quadrupole sources of various orientations 
and Di(u> 0 ) and Qij(u n ) are the dipole and quadrupole 
source strength respectively. In (2) k n is the wave number 
(jfc„ = Uo/aoo) and the source frequency u 0 is related to 
the observer frequency u through the familiar Doppler 
factor 

W = . TT~* — • (3) 


1 ~ M r COSd n 


When the quadrupole source is isotropic ( Qij = Q n for 
i,j = 1,2,3), it can be concluded that 


Qo ^ 4t tR a® (l — McCOsOoq) 5 dfjtoo 


(4) 


A Doppler factor correction for source volume effects 11 
has been utilized in (4). 


2.2 Three-Dimensional Geometrical Acoustics 


For an inviscid flow governed by linearized gas dy- 
namic equation, with mean velocity U(X) and sound 
speed a(X), the propagation of high frequency noise is 
governed by 


dpi 

ds 
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where Af, denotes the Mach component Ui/a , p, is the 
component of phase normal p and C(X) is the normal- 
ized sound speed defined earlier. Equations (5a) are nu- 
merically integrated along the ray subject to the initial 
conditions 

X = Xo at $ = 0 


„ ^ i M 

Po = AX + (AX-M--)^. (56) 

The ray speed A =| <fX | /ds is given as 

A = [/? 2 + (X - M) 2 ] - ’/C (5c) 


P 2 = 1 - | M | 2 

where unit vector X„ = (cos p, sin p cos 6, sin p sin 6) de- 
notes the direction of emission at the source (figure 2). 

Angle fi is measured with respect to a stationary 
frame of reference and is simply the vector sum of the di- 
rection of emission in the source coordinate and the mean 
velocity vector. For subsonic jets (| M n |< 1) this angle 
could clearly become as large as 7r, however, for a super- 
sonic flow (| M 0 |> 1) the upper limit on angle fi is less 
than 7t/2, resulting in a real value for the ray speed A 
in (5c). When the ray becomes straight A approaches one 
and a pair of angles ( 0 oo , <f>oo) are calculated for every pair 
of pre-selected {fi, 6). The far field mean square pressure 
directivity for a convecting multipole source of frequency 
u 0) accounting for source volume effects becomes 9,10 


~p2~ / P<*> q o 0 ,2, g4 ” 

M a ' 4 ttR ' p n sin 6 | | 


)x 


(1-^ Poo) 2 

(1-Jb Pn) 2 (l~^ P«) 




I (Po)iDi | 2 


( 6 ) 


P L = £ - )' i (p-Wi i 5 

1 *OQ P ° 

Here the source and observer frequencies are related as 
W = Wo /(l-^p„). (7) 

Uoo 

For an isotropic quadrupole source QijpiPj = Qo I P | 2 
and the sound field for a source of strength Q 0 {va) be- 
comes 
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( 8 ) 

In the far field Uoo is usually small. It is shown 
in the Appendix that A a /p„ is of the order of a~ a . It 
appears that equations (4) and (8) are essentially simi- 
lar except for the directivity factor d6/d<j > which is re- 
placed with its 3D counterpart in equation (8). In addi- 
tion, the 3D Doppler factor is influenced by the refraction 
phenomenon, as there is a unique phase normal p 0 at 
the source which results in propagation of a ray through 
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an observer point. It will be shown later that the sin- 
gularity in (3) in the direction of eddy Mach emission 
{Boo = cos -1 1 /Mr.) does not appear in (7). In fact we 
will compare the two expressions and assess the noncom- 
pactness constant of the modified Doppler factor. 

3. Numerical Results 

In order to integrate the propagation equations (1) 
and (5) it is necessary to specify the spatial dependen- 
cies of the mean jet velocity and temperature. A closed 
form formulation will naturally speed up the computa- 
tion of the gradients. However, when the jet is calculated 
numerically from the commonly available CFD codes, it 
is possible to integrate the CFD solution with an appro- 
priate interpolation routine (2D interpolation for axisym- 
metric jets and 3D for non-axisymmetric) and fit a muti- 
dimensional function on data at the grid points. Having 
found the tensor coefficients of the interpolating function, 
the directional derivatives at any point within the flow is 
readily available. 


3.1 Closed form Jet Model 


An axisymmetric self-similar Mach profile can be 
written in spherical coordinates as 


U(/?, 9) 
a(R, 9) 


fU . -(A)" 

(— )ci« 


( 9 ) 


where U is in the direction of 9, ct is the jet spreading 
angle and ( U/o) C l denotes the centerline value of the 
Mach number. For convenience, the velocity decay on the 
centerline is expressed as a function of the axial distance 
from the exit plane 10 . Further, in the case of a cold jet 
(constant stagnation temperature), the sound speed is ex- 
pressed using adiabatic flow relationship. For a subsonic 
jet with exit Mach number of U 3 / a ^ = 0.99, the radiated 
noise for an off-axis source was worked out in Ref. 10. In 
that exercise, the jet spreading angle was a = 10° and a 
3D ray-acoustic approach was employed to assess the di- 
rectivity of a convecting ring source at 7 diameters from 
the exit. 

For the purpose of comparison with the 2D predic- 
tions, we now let the jet spreading angle be very small 
(a = 1°), while the streamwise evolvement of the flow ex- 
pressed through n(Xi ) is kept intact ( see Ref. 10). Thus, 
in the context of 2D-GA, the flow is more or less a quasi- 
3D jet. This is due to fact that rays obtained according to 
(1) are representative of the flow characteristics on a span- 
wise plane at the source location. By moving the source 
to a new axial location the surrounding flow is allowed to 
change accordingly. 


3.1.1 Directivity Angles 


Figure 3 shows emission angle 6 v.s. azimuthal an- 
gle <poo with polar angle 6„ 0 indicated as a parameter 
for an off-axis source at x„/D = 7 and r„jD = 0.31. 
A proper selection of the target angle 9 no is essential 
in solving (1). There is a minimum allowable angle for 
any source point found by setting ql = 0, which gives 
B* = cos- 1 l/Co(l + Mo) as the boundary of the zone of 
silence. An attempt to solve the propagation equations 
with 6oo < 9* will trap the rays within the envelope of 
ql = 0 as depicted in figure 4. The boundary of the zone 
of silence for this source is at 53.7° and figure 4 was ob- 
tained for 9„o = 50.0°. When $«, > 9* but still close 
to 9*, the rays are confined to a narrow azimuthal an- 
gle (fino as seen in figure 3 for < 55° . This is due to 
a complete refraction of rays as will be shown clearly in 
3D-GA. In addition, it appears that within this range of 
Boo two rays could be arriving from a single source. How- 
ever, the 3D directivity angles indicate that the two rays 
indeed arrive at slightly different polar angles. Our 3D 
ray-acoustic computations for the same source point are 
shown in figure 5, where observer angles are given as a 
function of radiation angles at the source. Most of the 
radiated energy is limited to a narrow range of &oo , near 
the zone of silence (figure 5b). An enhanced version of 
figure 5b is given in 5c, showing slight variation in #oo 
with 6 resulting in a single ray arrival from each source 
point. Here the boundary of zone of silence was 52.7°, 
differing from the 2D predictions by only one degree due 
to non-parallel effects. In the limit, as the flow becomes 
parallel, 9oo becomes a function of /i only 


cos Boo — 


0,(1 - MS) 


( 10 ) 


and with n = 0 the boundary of zone of silence reduces to 
that of 2D ray-acoustic described earlier. For the purpose 
of completeness, we point out that figure 3 is, more or less, 
a direct conclusion of figures 5a and 5b. It is seen from 
figure 5b that if the variation of with S is neglected, 
parameter /i will be associated with a single 9 X (say 6m 
at 6 = 0). Eliminating cone angle /i will combine 5a and 
5b into figure 6 which shows a remarkable resemblance to 
the directivity angles of figure 3. 


3.1.2 Sound Pressure Level Directivity 


We continue our comparison of the 2D and 3D ray- 
acoustic methods with the computation of the noise di- 
rectivity for a convecting isotropic quadrupole source as 
given in equations (4) and (8). Although a prior knowl- 
edge of the source strength and its frequency is a re- 
quirement, for our purpose we limit the discussion to 
a point source or a ring source within an axisymmetric 
flow and leave out dependence on Q^u„. The convec- 
tion velocity {/„ in the most energetic part of the shear 
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layer is close to the jet mean velocity (Ref. 13) which 
leads one to believe that the high frequency noise scales 
as —3 power of the Doppler factor. In all the following 
computations, however, the convection Mach is expressed 
as a weighted average of the local and exit Mach num- 
bers (referenced with respect to ambient sound speed) as 
U^/aoo = O.SUa/aoo +0.311,-/(100. Figures 7a, 7b and 
7c show the sound pressure level directivity for an off- 
axis quadrupole point source at 7 diameters from the jet 
exit. On approaching the boundary of the zone of si- 
lence, the 2D directivity factor becomes infinitely large 
(figure 7a). A 3D ray-acoustic prediction (figure 7c), on 
the other hand, gives a finite value for the directivity fac- 
tor as both sin and | I become small. Further 

away from the zone of silence (0oo > 57°), there is close 
agreement between the two predictions. A different view 
of the spherical directivity for this source point is shown 
in figures 8a and 8b. 

The sound emission from a ring source is the direct 
result of integration of figure 7 for 0 < <j>oc < and is 
shown in figure 9. Here, the 2D calculation were termi- 
nated close to the zone of silence where the 3D predic- 
tions result in a sudden drop within the zone of silence. 
While the decay rate within the sound valley will decrease 
with increasing the jet spreading angle, still the high fre- 
quency approximation predicts a faster decay compared 
with measurements. Ref. 10 compares data with the 
ring source predictions for this jet at a spreading angle 
of a = 10°. 

3.2 Supersonic Jets 

The geometry considered for this exercise is an ax- 
isymmetric convergent-divergent nozzle. The flowpath is 
designed to generate a uniform and parallel flow at the 
nozzle exit for the design Mach number of 1.4 (figure 
10a). Both aerodynamic and acoustic data 14 were avail- 
able for comparison with predictions. For this hot jet 
the stagnation temperature was 1716 0 R which results in 
an exit Mach number of Mj = Uj/a <■» = 2.16. An ax- 
isymmetric version of the PARC code with Chien’s k-e 
turbulence model 16 was employed to solve the complete 
Reynolds-averaged Navier Stokes equations in conserva- 
tive law form. The predicted flow field 17 was shown to be 
in good agreement with data of Ref. 14. A third order 
spline fit was used to interpolate a surface on the CFD- 
generated flow field at the grid points (91 x60 points). The 
spline fit for the normalized sound speed C(X) is shown 
in figure 10b. For each one of the four variables of inter- 
est (sound speed and the three components of Mach vec- 
tor) the tensor coefficients of the interpolating functions 
were passed to a subroutine that computes the directional 
derivatives. 


3.2.1 Directivity Angles 

Figure 11 shows the 3D predictions for the boundary 
of the zone of silence for this supersonic jet as a function of 
source location. In a fine-grained turbulence, the acousti- 
cally active sources are associated with the kinetic energy 
of turbulence which is primarily distributed around the 
lip-line of the jet ( r 0 /D = 0.5). We will see shortly that 
for supersonic jets the radiation from any source point is, 
more or less, concentrated in a narrow range of polar an- 
gle near the zone of silence of that source. It is therefore 
reasonable to conclude that on the basis of figure 11 peak 
noise for this jet is in direction of 42° < 6 DO < 50°, in 
agreement with experimental data 14 . 

The radiation by turbulence convected at supersonic 
speed was termed by Phillips 18 as eddy Mach wave emis- 
sion. The Mach waves are believed to be emitted from 
the shear layer at an angle of 6 = cos' where /? is 

a constant generally taken as 0.7. This results in a Mach 
angle of 6 = 48.5° for the present jet. Another important 
feature of supersonic jets is related to source efficiency. As 
has been pointed out by Ffowcs Williams 19 the individual 
quadrupole sources essentially radiate as monopoles with 
no near field cancellation. This increased efficiency is ca- 
pable of converting between 0.1 to 1 % of the mechanical 
energy to noise as opposed to .01 % at angles beyond this 
range. 

The 3D directivity angle for an off-axis source at 
x 0 /D = 7 and r 0 /D — 0.5 is shown in figure 12. An en- 
hanced version of figure 12b is shown in 12c where the 
slight variation of Ono with <5 is noticeable. A distinct 
zone of silence at 6* = 48.3° is seen clearly. With equal 
increments of 5° in the emission angle /x, the radiated en- 
ergy is shown to be concentrated at a narrow angle near 
0*. As /i— + 0, there is a possibility for both and 
to become increasingly small. When this occurs the direc- 
tivity factor will approach § but it is possible to resort to 
higher order derivatives at the cost of an increase in the 
required precision in the directivity angles. 

To assess the performance of the 2D-GA, the directiv- 
ity angles for the same source point within the supersonic 
jet is shown in figure 12d. The zone of silence is given 
as 66.7°, not anywhere close to the previous prediction. 
Presumably the non-parallel effects play a major role in 
refraction of acoustic energy at supersonic speeds. 

3.2.2 Noise of Supersonic Jets 

In the last section, we briefly addressed the role of 
fine-grained turbulence in the noise of supersonic jets. It 
was shown that noise is primarily radiated in a narrow 
angle near the zone of silence. When the jet is not per- 
fectly expanded, shock-associated noise and screech tones 
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are generated due to interaction of turbulence with the 
shock cells. Here we concentrate on the mixing noise as 
described through equations (6) and (8). An exami nati on 
of equation (6) indicates that the monopole field Pjt, is 
of the order —3 power of the Doppler factor, as is equa- 
tion (8). Leaving out the source strength, both equations 
produce similar directivity patterns. 

Figures 13a & 13b show the spherical directivity of 
the radiated noise for an off-axis source, with parameter 
^ indicating the relative azimuthal angle between the 
source and the observer. As was pointed out earlier, a high 
clustering of rays occurs within a narrow polar angle at 
the boundary of the zone of silence, resulting in increased 
difficulty in the computation of the directivity factor. 

The role of refraction and its interaction with the 
Doppler factor for supersonic jets is demonstrated by com- 
paring the denominators of equations (3) and (7) in fig- 
ure 14a. Due to the singular behavior of the Doppler fac- 
tor, Ffowcs Williams 11 and Ribner 12 suggested a modified 
Doppler factor C"* = {(1 — M r cos #oo) 2 + a 2 M 2 } 0 5 . The 
noncompactness constant a is a semiempirical constant of 
the order of 0.3 — 0.55. Here we compute a directly from 

£^* = (1-— P„). (11) 

Goo 

The results are shown in figures 14b & 14c. It is seen from 
14b that a is angle dependent, approaching 0.4 for the 
most energetic part of the radiation field. A more accurate 
representation of <x is given in figure 14c as a function 
of the radiation angle p. The reader may recall that in 
our model the eddy convection Mach is a function of the 
source location. By moving closer to the axis ( r„/D = 
0.25), the Doppler factors compare as seen in figure 15a 
resulting in a = 0.33 shown in figures 15b k 15c. Overall, 
by investigating various source locations within the jet, it 
appears that the suggested value for a is indeed a good 
approximation. 

We show the directivity of the ring source in figure 
16a. The scaling of the directivity with power -3 of the 
Doppler factor is also shown in this figure as 



where the integrand is calculated at a point on the ring. 
Moving the ring source closer to the jet axis will result in 
directivity figure 16b. 

Our ring source predictions are also compared with 
data of Ref. 14 as seen in figure 17. It appears that the 
most important features of the noise directivity of this su- 
personic jet are captured through the ring source elements 
selected in this exercise. An integration of equation (8) for 


the entire jet can be worked out with proper knowledge 
of the source intensity spectra. 

4. Conclusions 

In this paper, applications of geometrical acoustics 
to propagation of high frequency jet noise was discussed. 
The propagation equations were solved using both 2D k 
3D ray-acoustic models. It was demonstrated that at sub- 
sonic Mach numbers, the 2D analysis results in satisfac- 
tory predictions outside the zone of silence, and therefore 
should be considered a viable substitute for the more com- 
plicated 3D computations. Unlike the subsonic jets, the 
parallel flow assumption proves to be inadequate for high 
speed jets. 

By computing the mean flow using a CFD code, it 
was shown that the directivity factor can indeed be cal- 
culated for jets of arbitrary geometries. Assuming that 
the sources of noise are embedded within the mean flow , 
we showed that refraction phenomenon modifies the su- 
personic Doppler factor significantly and results in a non- 
singular pattern in the direction of eddy Mach wave emis- 
sion. These arguments were used to assess the modified 
Doppler factor based on the line-of-sight method. 

The extension of the problem to non-axisymmetric 
jets seems relatively straightforward, although the in- 
crease in computational requirements suggests that per- 
haps a number of jet slices should be included in the vol- 
ume integration. In addition, such an attempt will also 
require the magnitude of the source strength of the non- 
axisymmetric jet. 
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Appendix 

We approximate X a assuming that the mean flow is 
in X\ direction, i.e., M = Mi, and that the relative az- 
imuthal separation between the direction of emission and 
the observer is zero, i.e., X„ = ( cos p , sinp , 0). From 
equation (5c) we have 


A7 l = — x/I (,41) 

<*oo 

The term M„ sin /i is usually smaller than one. Even when 
the jet becomes supersonic , i.e. M„ > 1, most of the 
radiated energy is limited to small values of the emission 
angle p; hence A7 1 ~ Using the one-dimensional gas 

dynamic relation for an adiabatic gas we have 

(.42) 

Poo ^00 

With 7 = 1.4 it is seen that (^ )A 8 ~ (^) 8 . 
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Figure 6.—Combinatk>n of figures (5b) and (5c) when variation 
of 0. with 6 is neglected. 
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Figure 9. — Noise directivity for a ring source at X</D = 7 and rj/D = 0.31 . 
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Figure 12. — Directivity angles of the supersonic jet for an off-axis source at X{/D = 7 and rg/D = 0.50. Figures {a) and (b) 
are the 3D-GA predictions, (c) is the enlarged (b). (d) is the 2D-GA prediction. 
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Figure 1 4. — (a) Doppler factor of the supersonic jet for the source of 
figure 13. (b) Empirical constant of the modified Doppler factor for 
the source of figure 1 4a. (c) Em pineal constant of the modified 
Doppler factor for the source of figure 1 4a 



Figure 15.— (a) Doppler factor of the supersonic jet with the 
source point at XyD = 7 and r</D = 0.25. (b) Empirical 
constant of the modified Doppler factor for the source of 
figure 15a (c) Empirical constant of the modified Doppler 
factor for the source of figure 1 5a 
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Figure 1 6. — Noise directivity for a ring source within the supersonic jet at 
X(/D = 7. (a) r(/D = 0.50. (b) r^D = 0.25. 
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Figure 1 7.— Comparison between data and high frequency approximation. Predic- 
tions are for a ring source at X</D = 7 . 
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